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Abstract 



We analyse derivative securities whose value is not a deterministic function of an under- 
lying which means presence of a basis risk at any time. The key object of our analysis is 
conditional probability distribution at a given underlying value and moment of time. 

We consider time evolution of this probability distribution for an arbitrary hedging strat- 
egy (dynamically changing position in the underlying asset). We assume log-brownian 
walk of the underlying and use convolution formula to relate conditional probability distri- 
bution at any two successive time moments. It leads to the simple PDE on the probability 
measure parametrized by a hedging strategy. For delta-like distributions and risk-neutral 
hedging this equation reduces to the Black-Scholes one. We further analyse the PDE and 
derive formulae for hedging strategies targeting various objectives, such as minimizing 
variance or optimizing quantile position. 
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1 Introduction 



In the classical theory of derivatives a number of restrictions is conventionally imposed. 
Risk free rate of interest, growth rates and volatilities of log- normal processes are assumed 
to be constant, short selling of securities is permitted, transaction costs, taxes and div- 
idents are absent, securities are perfectly divisible, trading is continuous, and there are 
no arbitrage opportunities [1] . Attempts to overcome these restrictions may involve some 
basis risk. In competitive markets the price is such that a participant may receive a profit 
or suffer a loss. We expect that market participants are prepared to price basis risk in 
a form of profit/loss probability distribution function (PDF) according to their portfolio 
needs. Our problem here is to develop a technique for computing PDFs of derivatives for 
any hedging strategy. 

Risk-neutral valuation is a general approach to studying derivatives [1]. It basically says, 
that in order to compute a present value of a contract (European, for example), one 
must average the cash flow function at maturity over equivalent risk-neutral martingale 
measure Ai[S(t)] on a space of random walks of the underlying asset [S(t)\ and make an 
appropriate discounting 

F(S,t) = e r ^E M [F(.,T)} (1) 

This prescription follows from the no- arbitrage argument. The effective risk- neutral mea- 
sure M comes as a result of gauging of an objective measure M[S(t)] by the riskless 
hedging strategy (f> (S, t) 

(M, O ) — ► M (2) 

Riskless hedging is merely a replication process used for correct pricing. In addition 
to replicating riskless portfolio an investor may choose to explore innovative hedging 
strategies. Then investors should be provided by a tool which would enable them to 
choose hedging strategies serving best their personal objectives. 

So, what if a generic hedging strategy ^ 0o is chosen? This would lead to an effective 
measure M.^ which is generally neither risk-neutral nor martingale. We shall solve the 
following problem. 

Given the objective measure Ai (say, log-normal with drift ji and volatility a) 
and a hedging function </>, derive an effective measure M.^. 

It is worthwhile noting that this problem is also related to pricing in incomplete markets 
[2-4], when O is not unique. 

As long as effective measure Ai<p is not riskless a value of the contract F is clearly a random 
number. Therefore, we should describe M^, by a conditional probability distribution 
function (PDF) V^(F\S,t) for a given S (underlying value) and t (time): 

/oo 
dF V+(F\S,t) = 1 (3) 
-oo 
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Then equation (1) takes more general form 



V*(F\S,t) = e riT - l) E M ^[V(Fe r ^\ ■ ,T)} (4) 

Time evolution of V^(F\S, t) can be described by the PDE which we derive in the next 
section. It can be viewed as generalized Black-Scholes equation. In fact, we show that it 
coincides with the Black-Scholes equation in the risk free limit. In the sequel we omit the 
(f> index. However, we always imply that a PDF depends on 0. 

2 The PDE on a PDF 

2.1 One Time Step 

Let us consider the time interval (t, t'). Assume that the underlying stock value at the 
beginnig S(t) = S and conditional PDF V(F\ ■ ,t') at the end of the interval as well 
as a risk neutral discount rate r are known. Markovian transfer matrix p(S' , t'\S, t) will 
be considered log-normal in the sequel. We should also assume that hedging position <fi 
is unchanged during the interval. We consider a position where the contract is sold and 
shares are purchased. Portfolio consisting of shares and one short contract has a present 
value X = F — <pS. A PDF of the portfolio value X at the moment t is simply expressed 
by the following convolution 

V(X\t) = J dS' p(S',t'\S,t)V(X + (f)S'\ S',t') (5) 

Now suppose that a value of the potfolio at the moment t' is X with some probability 
p. Then clearly a contract value at the moment t is F = X e r ('~'') + <f>S with the same 
probability. This simple reasoning leads to the useful relation between the contract value 
statistics at the moments t and t' 

V(F\S, t) = e rit '- l) J dS' p(S', t'\S, t) V((F - (/>S)e r ^ + </>S'\ S', t') (6) 

So we did a remarkable loop in time when deriving (6). First, we used forward transfer 
matrix on S to mix S and F statistics at t' and obtain statistics of X, Eq(5). Then 
we went back in time discounting portfolio statistics and subtracting stock part of the 
portfolio at t. 

Another important lesson which we learn from this exercise it that statistics of a contract 
value is closely related to a hedging strategy. To find out how valuable a contract is one 
needs to examine it's statistical behaviour under optimal hedging. We should return to 
this point in the next section. 
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2.2 Continuous Limit 



We now take a continuous limit assuming that the time step (t, t' — t + dt) is small and 
therefore, the kernel 

is sharply peaked as compared to the typical support of V(F\S, t). Then expanding r.h.s. 
of (6) to the first order in dt one finds [5] 

f + r A «F - + „ S + % ) + ^ + 2^ + §*) = (8) 

Second term in the l.h.s. corresponds to continuous discounting of the portfolio, third 
term is responsible for the drift of the underlying and fourth one is a diffusion term a 
la Fokker-Plank. It is instructive to derive equations for mean, conditional variance and 
skewness of the contract value. Let us define 

F(S, t) = E V [F] , V(S, t) = E V [F 2 } - F 2 , Q(S, t) = E P [F 3 ] - 3E V [F 2 }F + 2F 3 . (9) 

Differentiating (9) with respect to time and using (8) gives 
f)F clF 1 f) 2 F 

m + » s os + r 2s2 8<p- rF = ^- r ^ ^ 



dV dV 1 2 2 d 2 V 

-m + » S ds + 2 aS ds^- 2rV = 



« 2 <i 2 ( dP 



and 



dQ JO 1 2 2 d 2 Q 



(12) 



Pre-factor 2 (3) in the term 2rV (3rQ) accounts for the fact that the variance (skewness) 
is measured in squared (cubed) units of currency. Equation (10) reduces to the Black- 
Scholes equation on the average F if one chooses the hedge <fi = OF /OS. If one adds to 
it final condition on variance V(S,T) = then (11) gives V(S, t) = 0, i.e. randomness 
disappears at any moment t. Thus we obtain the classical European contract dynamics 
special case. 

PDE's for the momenta E-p [F n ] can be easily derived in the same fashion. Lower momenta 
E-p[F n ~ l ] and E-p[F n ~ 2 } will enter the nth equation as long as we have got first- and 
second-order partial derivatives in (8). 
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2.3 Path Integral Solution and Effective Measure 



Deriving equation (8) we had in mind that a PDF is a smooth function of it's arguments 
(otherwise partial derivatives would not be well defined). Unfortunately, it is not always 
true. The PDE is not applicable if we have to deal with binary, delta-like or any other 
singular distribution at maturity. However, it is possible to derive closed form integral 
solution (evolution kernel) which would make sense for any reasonable final conditions 
V(F\S,T). 

Suppose that we are interested in a PDF at initial moment t. Let us divide the time 
interval (t,T) in N little segments (tk,tk+i) such that to — t and = T. We should 
move back in time - from maturity to the present moment - applying backward transfer 
matrix (6) at each step. Introducing notations S(tk) = Sk and (f>(tk) = (fik one can rewrite 
(6) as 

V(F\S k ^,t k ^) = JdSkT^S^lSkMFlS^h). (13) 

The T-operator acts on both arguments of V - on S as a p-matrix and on F as a shift. 
After N successive backward steps we get 

V(F\S,t) = J dS 1 dS 2 ...dS N T (P (S\S 1 )T (P (S 1 \S 2 )...T (P (S N - 1 \S N )V(F\S N ,T). (14) 

In the N — > oo limit the p-matrices accumulate to the log-brownian measure 

J DM\ S (t)=s = J dS 1 dS 2 ■ ■ ■ dS N p(S,t\S 1 ,t 1 )p(S 1 ,t 1 \S 2 ,t 2 ) ■ ■ . p(£jv-i, *n-i|-SV, *jv)> 

(15) 

whereas stepwise shifts of F produce an integral shift 

V[S,<j>] = e~ rT £ <P(dS(u) - rSdu)e ru (16) 

defined by the path S(u) and the hedging function (f>(S,u). So, we get a compact answer 
for the PDF 

V(F\S,t) = e riT -^ J DM[S] V{e r{T ~ l \F + ^[S,</>]) \ ■ ,T) (17) 

It's a right time now to recall the formula (4). Effective measure M.^ on the space of 
random walks [S] turns out to be a differential operator acting on V 

J DM^[S] = J DM[S] exp{^[S,<t>] d F ). (18) 

It should be stressed that formula (17) does nor require any special properties of the 
probability measure like smoothness or continuity. It perfectly works for any probability 
distribution. As an example let us consider a European call option with random strike. 
Say, the strikes are K\ and K 2 with probabilities p and 1—p respectively (hybrids or dual 



5 



triggers are contracts of this type). Maturity distribution for such an option is a linear 
combination of two 5-functions and equation (17) immediately gives 

F*(F|S,0 = e r{T - l) J DM[S] {p S^-'^F + *[S, 0]) - (S(T) - K ± ) + ] + 

(1 - p) Sie^iF + 0]) - (S(T) - K 2 ) + }}. (19) 

3 Optimal Hedges 

In this section we show how to optimize different characteristics of a PDF. First, let us 
consider variance at the moment t: V(S,t). Suppose that the goal of a hedger is to 
optimize this variance by using the best strategy [6]. One should minimze variance as a 
functional of 0(S",t'), where t' G (t,T). A convenient Green function 

PA(S',t'\S,t) = e Ar ^ p(S',t'\S,t) (20) 

turns out to be useful to express evolution equation solutions whose dimensionality with 
respect to a numeraire is A. For example, solution of the equation for variance (11) can 
be presented as 

V*(S,t) = J dS'dt' p 2 (S',t'\S,t){v(S / ,t , )5(t'-T)+a 2 S' 2 ((f ) (S',t')-ds'F(S',t')) 2 } (21) 

Variance minimizing hedge 0* solves the variational equation SV^/S4> = 0, which imme- 
diately gives <f>*(S,t) = dsF(S,t). Eq. (10) for the variance minimizing hedging strategy 
0* is nothing but the Black-Scholes equation on the mean. 

Final condition to be used here is F(S,T) = Ep[F]\ t= T- Thus variance minimizing hedge 
is given by 

,*,q +\ 9F BS (S : t) 

<? [S,t) = . (23) 

Minimal variance V^~* can be found then as a solution of the homogenous PDE 

Another important characteristic of the statistical distribution which we consider is a 
quantile position [7] F q (S,t) defined by 

/ dFV{F\S,t)=q (25) 

J — oo 
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Optimal hedging 0** now is such that max^\F^\ = F^** . Similar variational equation 
5 F^ /5(f) = can be solved for this case as well. It leads to 

d F lnV* \ F =F q (s,t) 

where V^** is the solution of (8) for = 0**. Although the system of two equations 
(8) and (26) looks nonlocal in time, it makes perfect sense once we adopt the stepwise 
approach. Namely, we assume again that the time is discrete (t,ti,t2, ■ ■ - t^ — T). Then 
knowing V(F\S, tk) and therefore F q (S,tk) one can compute (f)**(S,tk-i) using formula 
(26) and V(F\S, tk-i) using (13) (discrete version of (8)), etc. What we describe here is a 
dynamical programming procedure which would result in a hedging function maximizing 
a quantile position. Such a procedure will be implemented numerically and presented in 
a separate paper. 



4 Conclusion 



In this paper we have indicated an approach which allows one to study financial derivatives 
depending on the method of their hedging. We think that this approach can be useful for 
developing optimal trading strategies in the framework of portfolio management. This 
approach can also be helpful for quantitative analysis of contracts which value at maturity 
is random rather than deterministic. Randomness may come from various sources such 
as default, contingency on non-traded indices, or uncertainty of the statistics of trade- 
able underlyings. This method also leads to analysis of incomplete markets, and helps to 
develop preferences regarding investment strategies (the case V(F\-,T) = 6(F)). In the 
latter case, an intersting problem is to find the hedge which replicates a preferred proba- 
bility distribution. One way to do this is optimize the Kullback-Leibler distance between 
the replicated and target distribution functions [8]. The work on application of objective 
measure analysis to American options and fixed-income instruments is in progress. 



5 Acknowledgment 

S.E. is grateful to Jay Blumenstein for discussions, to Chunli Hou for help with the 
literature on mean- variance hedging, and to Dajiang Guo for applying these ideas to 
pricing derivatives on S&P 500 in objective measure. I. V. thanks Sasha Migdal for support 
and sharing his mathematical ideas. We are thankful to Marco Avellaneda for his valuable 
suggestions and encouragement. 



7 



References 



[1] J. Hull, Options, Futures and Other Derivative Securities, Prentice Hall, Englewood 
Cliffs, NJ, 1993 

[2] H. Follmer and D. Sondermann, Hedging of non-redundant contingent claims. In Con- 
tributions to Mathematical Economics (W. Hildenbrandt and A. Mas-Colell, eds.), 
1986, 205-223 

[3] D. Duffie and H.R. Richardson, Mean-variance hedging in continuous time, Ann. 
Appl. Probab. 1, 1991, 1-15 

[4] H. Follmer and M. Schweizer, Hedging of contingent claims under incomplete infor- 
mation. In Applied Stochastic Analysis (M.H.A. Davis and R.J. Elliot, eds.), Gordon 
and Breach, London, 1991, 389-414 

[5] S.E. Esipov, unpublished 

[6] M. Schweizer, Mean-variance hedging for general claims, Ann. Appl. Probab. 2, 1992, 
171-179 

[7] H. Follmer and P. Leukert, Quantile Hedging, preprint, 1998 

[8] M. Avellaneda, C. Friedman, R. Holmes, and D. Samperi, Calibrating Volatility 
Surfaces via Relative-Entropy Minimization, preprint, 1997 



8 



